Generation of maximally entangled charge-qubit arrays via a cavity mode 
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We describe an assembly of TV Cooper-pair boxes (CPB) contained in a single mode cavity. In 
the dispersive regime, the correlation between the cavity field and each Cooper-pair box results 
in an effective interaction between CPBs that can be used to dynamically generate maximally 
entangled states. With only collective manipulations, we show how to create maximally entangled 
quantum states and how to use these states to reach the Heisenberg limit in the determination of a 
spectroscopy frequency. This scheme can be applied to other types of superconducting qubits. 

PACS numbers: 03.67.Mn, 03.67.Lx, 85.25.Hv 



The description of the interaction between atoms and 
quantized modes of the electromagnetic field in a cavity is 
called cavity quantum electrodynamics (cQED). The first 
experimental studies used flying Rydberg atoms in an rf 
resonator. With the advent of quantum computing sev- 
eral other implementations were developed to mimic the 
quantum properties of atoms. Among those, solid-state 
implementations are especially interesting since they offer 
several advantages over real atoms: these artificial atoms 
properties can be tailored and their number and location 
is fixed. We concentrate our discussion on a supercon- 
ducting circuit called a Cooper-pair box (CPB). Several 
experiments have shown that these circuits behaved like 
quantum two-level systems 0, 0, 0, H 01 • An architecture 
has been proposed to realize an on-chip cQED experi- 
ment using a CPB as the artificial atom strongly coupled 
to a one-dimensional cavity Q. This architecture was 
implemented and the strength of the coupling was shown 
to be indeed stronger than the the different decays con- 
stants so that the vacuum Rabi splitting was observed 
0- Subsequently the ac-Stark shift was measured us- 
ing a quantum non-demolition technique (QND) 0] and 
the decoherence time T 2 evaluated from Ramsey fringe 
experiments 0. 

The possibility of entangling two Cooper-pair boxes 
in a cavity QED scheme was described in reference 0. 
Beside its fundamental interest and applications to quan- 
tum information processing, entanglement offers the ad- 
ditional advantage of allowing improved sensitivity in 
a quantum-limited measurement. Interferences between 
two different polarizations (modes) of three and four pho- 
ton entangled states have been observed 0, . Spec- 
troscopy performed on an assembly of three beryllium 
ions also demonstrated a similar improvement in fre- 
quency estimation ^3] • Using entanglement in solid-state 
implementations to enable such Heisenberg limited mea- 



surements could revolutionize sensor technology with, for 
instance, electrometers and magnetometers. 

We propose to use the photon-qubit interaction to cre- 
ate an effective interaction between all the CPBs. We 
show how this interaction can be used to generate maxi- 
mally entangled states and in turn how to use these states 
to beat the standard quantum limit when measuring the 
Ramsey frequency. 

The Hamiltonian of a single Cooper-pair box can be 
expressed as [l3j : 



Hr 
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(1) 



The <7fc are the usual Pauli matrices (k = x, y, z) ex- 
pressed in the charge basis {|0), |1)}. In the linear ap- 
proximation, the diagonal term, or bias, depends lin- 
early on the gate charge n g , B z (n g ) — 8E c (l/2 — n g ) 
where E c is the charging energy. So adding an increment 
Sfig to the gate charge n g results in a new Hamiltonian 
HQ(n g +6hg) = HQ(n g ) — (8E c Sh g /2)a z . We assume this 
approximation to be valid in a limited range < n g < 1 
(in units of 2e). The tunneling matrix element B x is the 
Josephson coupling Ej. In reference to the case where 
this quantity is zero (and therefore the Hamiltonian diag- 
onal), the point n g = 1/2 is called the degeneracy point. 
The linear approximation is best justified in the vicinity 
of this point. For charge qubits, E c > Ej, so that the 
Hamiltonian is approximately diagonal in the charge ba- 
sis at n g = and n g = 1. In the remainder of this article, 
we work at the degeneracy point. 

We now consider the case of a Cooper-pair box placed 
in a single mode cavity and subject to a quantized field 
Sh g = n c (a^ + a) on its gate. We are interested in the 
strong coupling limit, g > k, and therefore we neglect the 
cavity decay rate n (and the decay of the qubits excited 
state). We express the Hamiltonian in the eigenbasis of 
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and use the rotating wave approximation (RWA): 

H = cj c (a f a + 1/2) + P./y + |(a f o-- + acr+) (2) 

Here, we define the coupling constant g = —AE c n c . 
We can use the above Hamiltonian to describe several 
Cooper-pair boxes present in the cavity when they are 
not directly coupled. Assuming each CPB experiences 
the same coupling to the cavity field, we can rewrite: 

H N = uj c {o) a + 1/2) + EjS z + g{a) S- + aS+) (3) 

We used the notation S z = | J2iLi a zi > N being the num- 
ber of qubits (N > 1). Even though the above Hamilto- 
nian is in principle valid for any N, we work far away 
from the thermodynamic limit N — > oo, and assume 
N < u! c Ej/g 2 . Equation J2J is also known as the iV-atom 
Jaynes-Cummings Hamiltonian. We use the eigenstates 
of S z to label the total quantum states \M) Z (Dicke states 
\J,M) Z with J = N/2). 

In the case of a single atom, the exact diagonalization 
can be performed in the so-called dressed-states basis [bH ] 
and leads to the energy spectrum E±^ = [n + l)w c ± 
\ y/ig 2 (n + 1) + A 2 . The detuning A is defined as Ej - 
uj c . In the large detuning regime g/A -C 1 the Af-atom 
Jaynes-Cumming Hamiltonian J2J can be approximately 
diagonalized with the transformation : 

U = exp(g/A(aS+-a^S-)) (4) 

We calculate H' N = UHnU^ and get 




We define x = .9 2 /^- The shift in the resonator frequency 
(linear in photon number h = a^a) is known as the ac- 
Stark shift. In the second term, the shift in the atomic 
frequency is known as the Lamb shift. The last term 
shows that a correlation x$ z exists between qubits when 
N > 1. 

Hamiltonian J5J contains the nonlinear interaction 
H S z — X^ 2 - It is well known that this so-called one-axis 
twisting term can produces a spin-squeezed state when 
applied to a coherent spin state ^f|. A coherent spin 
state is an eigenstate of the angular momentum operator 
S n = n ■ S with an eigenvalue N/2 [n is a unit vector 
pointing in the direction (#,</>)]. The time evolution of 
such an interaction term applied to an initial state per- 
pendicular to the z axis, for a time t sz — n/2\, produces 
a maximally entangled state (assuming N even. If N 
is odd, another rotation should be added in the process 

My 

|*!> = exp(-iH sz t sz )\g) = j=(\g) +* 1 + JV |e)) (6) 



we used \g) = e~ %s ^/ 2 \ - N/2) z and |e) = e~ lS ^\g). 
Thus, Hamiltonian (J5J can be used to dynamically pro- 
duce maximally entangled states. 

Provided the rotation can be made in a time substan- 
tially smaller than the entangling time, we propose to 
start from the ground state | — N/2) z , perform a n/2 ro- 
tation, say around the axis x, then let the system evolve 
freely during t sz = ir/2x and finally perform another tt/2 
rotation: 

U N = e +iS ^/ 2 U sz (t sz )e-^/ 2 (7) 

Un can be seen as a generalized tt/2 Ramsey pulse. Start- 
ing from an initial state | — N/2) z , it produces a superpo- 
sition state -^(|-A72) 2 + |+A72) Z ). A method has been 

proposed recently |l2l | in the context of ion traps, where 
one collective measurement demonstrates the use of max- 
imally entangled states to reach Heisenberg-limited spec- 
troscopy. It generalizes the idea of the Ramsey fringe 
experiment to particles. First a n/2 pulse, U^, is ap- 
plied to the system. Then the system acquires a phase 
during a free evolution period that lasts until another 
tt/2 Ramsey pulse is applied: 

\* final) = U N e~ lS ^U N \ - N/2), (8) 
= -i sin(^)| - N/2) z + cos(^)| + N/2) z (9) 

The probability for the system to collapse in the state 
| + N/2) z = | t) is then simply P T = |(1 + cos{Nip)). 
The phase ip can be tuned by controlling the free evo- 
lution time T in between the entangling sequences Un- 
Alternatively, the phase can be controlled by adjusting 
an offset phase in the n/2 pulses of the last Un sequence 
(l2^ . From the expression of Pf, we see that it is pos- 
sible to evaluate the phase with an uncertainty, at best, 
equal to 1/N. This so-called Heisenberg limit represents 
an improvement of l/y/N over the standard quantum 
limit (obtained for a assembly of N atoms randomly cor- 
related). This method has been proposed to improve 
atomic frequency measurements |l7|. For Cooper-pair 
boxes, the Josephson energy is analogous to the atomic 
frequency. So the method proposed here uses entangle- 
ment to improve the accuracy on the measurement of Ej 
with spectroscopy. 

The vast majority of experiments performed with a 
CPB uses a SQUID geometry so that the Josephson en- 
ergy Ej max cos(7T$/$o) can be tuned with an external 
magnetic flux $ ($o is the flux quantum). Hence, an 
improvement in Josephson energy determination results 
in an improved flux measurement (if the energy Ej ma7 . 
independently measured). Since Cooper-pair boxes are 
charge qubits, these devices are not effective magnetic 
flux sensors (electric charge and magnetic flux are two 
quantum conjugate variables). Therefore, we emphasize 
that the improvement in the flux determination is only 
relative (to the case where CPBs wouldn't be correlated). 
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Instead of the fluorescence method usually used in ions 
traps to read the probability Pf (or PjJ, the QND mea- 
surement proposed and used in reference H, S S is 
more appropriate to our scheme. Beside the term xSj, 
equation (JSJ contains another term, 2xS z a { a 1 describing 
the interaction field-qubits. Because S z commutes with 
this interaction and the Hamiltonian J5j, a quantum non- 
demolition measurement can be performed and S z mea- 
sured. If a continuous wave signal at the resonator fre- 
quency lu c is applied on the system, the term correspond- 
ing to the ac-Stark effect will introduce an equal but op- 
posite phase shift for the two components |±iV/2) z . If we 
start from an initial state (a| + N/2) z N/2) Z )\Q) 
where \6) is a coherent state, the time evolution corre- 
sponding to the interaction (g 2 /A) a^a a z will create a 
state a\9,+N/2) z +j3\ -9, -N/2) z where tanfl = 2. 9 2 /kA 
@ (k is the cavity decay rate u) c /Q). The phase differ- 
ence between the two states | ± N/2) z can be measured 
by a transmission measurement of the photons through 
the cavity. The probability Pj is extracted from the time 
dependence of 9. This measurement can be either contin- 
uous (weak) or pulsed (stronger). The scheme proposed 
here is scalable to any N (as long as the qubits are only 
coupled together though the cavity mode). 

For a 24 mm long cavity Q , we can reasonably assume 
that an assembly of four CPBs spaced by 100 /zm from 
each other, will experience a homogeneous field at the 
central antinode and at the same time that the CPBs 
will not be directly coupled. For a minimum detuning 
of approximately 100 MHz and a coupling strength of 17 
MHz, the time needed to generate a maximally entan- 
gled state would be around 346 ns. The squeezing time 
obviously has to be smaller than the decoherence time. 
The time T2 has been evaluated to be 500 ns in the case 
of one qubit @ . The squeezing time could be decreased 
further by increasing the couplings strength so that the 
condition of large detuning j C A for which eq. ijjjj is 
valid would still hold. In order to minimize the effect of 
the measurement, the probability Pj can be read after 
the last 7r/2 pulse @. We conclude that with existing 
technology it is possible to create maximally entangled 
states between N (four) charge qubits in a superconduct- 
ing cavity QED and improve by a factor ^/N (y/2) the 
accuracy on the measurement of the Ramsey frequency. 

The method we proposed can be simplified further. 
The general sequence iJBJ requires four ir/2 pulses (two 
per general Ramsey pulse Un)- This number can be de- 
creased by half. For this purpose, we examine the situ- 
ation of atomic type spectroscopy |l4| performed on the 
system at the degeneracy point. The effect of a classical 
periodic drive on the gate can be described by adding 
e(e~ tuJt S + + e^'5_) to equation ©. In the rotating 
frame at to: 

H u = 5\S Z + 5 2 a) a + g(a^S- + aS + ) + 

e(S + + S.) (10) 



with S\ = Ej — lj and 82 = lo c — lj. When the drive is 
resonant with the qubits, 8% = 0, this interaction term 
can be expressed in the interaction picture : 

H[ = g ( Sx + S '~ iS v e *2. t _ S z + iS v e _ j2£t)atgj52t 
+h.c. (11) 

We restrict ourselves to the weak drive limit 82 > £■ Us- 
ing a method described by D.F. James 0, we average 
this interaction over a time T 3> ir/e and derive the ef- 
fective Hamiltonian: 




+ (f) 2 ^2i " S * ] + + 2 ^ a)S ^ 

+ (f Ts^TYe [(52 - sl) (1 + 2atfl)&] (12) 

The advantage of this expression over equation JSJ is that 
it is expressed in a basis orthogonal to the original z axis 
so that the one axis twisting term is now proportional 
to S 2 . Starting from the ground state | — N/2) z and 
turning on the interaction H c g for a time approximately 
equal to (7^2) / '(4g 2 ) will produce a maximally entangled 
state. So, with the same initial state, only one manip- 
ulation is needed to generate a GHZ state with a clas- 
sical periodic drive whereas two n/2 pulses are needed 
with the expression The coupling g should also be 
smaller than the detuning 82 so that the entangling time 
is long enough compare to the average time T after which 
(|12fl is valid. This effective Hamiltonian is valid for any 
system described by a N-atom Jaynes-Cummings Hamil- 
tonian subject to a weak period drive resonant with the 
qubits. Shifts analogous to one found in the Hamiltonian 
(JSJ , in the atomic and especially the resonator frequency 
are present in 1)12(1 . Thus, the same kind of quantum 
non-demolition measurement as proposed in [6j can be 
performed. 

In this work, we described a collection of N Cooper- 
pair boxes contained in a single mode cavity. At the de- 
generacy point, we showed that the diagonal Hamiltonian 
contains an effective qubit-qubit interaction mediated by 
the cavity mode, of the form x^ z - This interaction can be 
used to prepare maximally entangled states. We adapt 
a method used in ion traps to demonstrate the use of 
maximally entangled state to reach the Heisenberg limit 
in the Ramsey frequency determination. We estimate 
that such an experiment is possible with existing technol- 
ogy. Additionally, we showed that the iV-atom Jaynes- 
Cummning Hamiltonian, with a weak classical drive on 
resonance with the qubits, exhibits the same properties 
as the same Hamitonian without a drive, but expressed in 
an basis orthogonal to the eigenbasis of the system with- 
out drive. This approach decreases the number of pulses 
needed in our scheme to demonstrate entanglement. The 
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results presented in this work are valid for two-level sys- 
tems that can be described with the Hamiltonian Q . For 
instance, it can be applied to an assembly of rf-SQUIDs. 
This scheme is especially interesting for solid state imple- 
mentation since it only requires collective manipulations. 
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